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Dirac and Gor'kov spectra in two color QCD with chemical potential 

Elmar Bittner a , Maria-Paola Lombardo b , Harald Markum a , and Rainer Pullirsch a 
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We analyze the eigenvalue spectrum of the staggered Dirac matrix in two-color QCD at nonzero baryon density 
when the eigenvalues become complex. The quasi-zero modes and their role for chiral symmetry breaking and the 
deconfinement transition are examined. The bulk of the spectrum and its relation to quantum chaos is considered. 
A comparison with predictions from random matrix theory is presented. An analogous analysis is performed for 
the spectrum of the Gor'kov representation of the fermionic action. 
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1. Chiral and diquark condensate 

The eigenvalues of the Dirac operator are of 
great interest for important features of QCD. The 
Banks-Casher formula jy relates the Dirac eigen- 
value density p(X) at A = to the chiral con- 
densate, £ = \(ipip)\ = lim E _>o hmy^oo irp(e)/V. 
The microscopic spectral density, p s (z) = 
limy^oo p (z/V"E) /VS, should be given by the 
appropriate result of random matrix theory 
(RMT) 0, which also generates the Leutwyler- 
Smilga sum rules || . 

A formulation of the QCD Dirac operator at 
chemical potential p ^ on the lattice in the 
staggered scheme is given by 



M xv = — 
' y 2a 



y"l [Uv(x)r) v (x)8y !3>H , - h.t 



U { (x)r, t (x)e^ v ^-Uj(y) m (y)e-^ X J \ , (1) 



with the link variables U and the staggered phases 
rj. The bilinear form of the fermionic action 
Sferm = ipMif> permits only to calculate the one- 
point quark-antiquark function {ipip). A possible 
extension of the Hilbert space is provided by the 
so-called Gor'kov representation of the action 



Sferm = C0> 



JT2 
\M tr 



j T 2 



(2) 



with external source j , which allows to extract the 
diquark one-point function (ipip) B. 



The expression for the diquark condensate 
given in [Q can be reformulated as < ifjij) > = 
lim^ limy^ oo (l/(2F)E„(Af + j)" 1 ), where 
the are eigenvalues of the matrix 



K 



-ir 2 M 4 



(3) 



2'^" 3 1 
whose dependence on the diquark source j is via 
jl. The \^ then play, in the study of diquark 
condensation, the same role as the eigenvalues of 
the Dirac operator in the case of chiral conden- 
sation. Finally we note that K is related to the 
"Gor'kov operator" in Eq. (Q) by a multiplication 
with a constant matrix. In the following we will 
concentrate on the analysis of the Gor'kov spec- 
trum. 

Our computations with gauge group SU(2) on 
a 6 4 lattice at f3 = A/g 2 = 1.3 have N f = 2 flavors 
of staggered fermions of mass m = 0.07. For this 
system the fermion determinant is real and lat- 
tice simulations of the full theory with chemical 
potential become feasible. They exhibit a phase 
transition at p, c k, m,/2 « 0.3 where the chiral 
condensate (nearly) vanishes and a diquark con- 
densate develops ||. 

In the upper row of plots in Fig. [l] we present 
the behavior of the eigenvalues of the Dirac opera- 
tor with increasing chemical potential. A dilution 
of the eigenvalues around the origin is observed. 
The lower row of plots depicts the eigenvalues of 
the "Gor'kov operator" for j = as a function 
of p. One clearly sees that their density at the 
origin increases except for large p. 
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Figure 1. Upper plots: Complex eigenvalues of the Dirac operator at various values of fi for a typical 
equilibrium configuration of two-color QCD (horizontal axes = real parts, vertical axes = imaginary 
parts, in units of 1 /a). Lower plots: Same for the "Gor'kov operator". 
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Figure 2. Upper plots: Density p(y) of small eigenvalues of the Dirac operator for two-color QCD on a 
6 4 lattice from fj, = to 0.4 (left). The loss of quasi-zero modes is accompanied by a vanishing of the 
chiral condensate. Chiral condensate extracted by different methods (right). Lower plots: Similar for the 
"Gor'kov operator" (see text). 
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In the upper left plot of Fig. || we compare the 
densities of the small eigenvalues of the Dirac op- 
erator at p = to 0.4 on our 6 4 lattice, averaged 
over at least 160 configurations. Since the eigen- 
values move into the complex plane for p > 0, a 
band of width e = 0.015 parallel to the imag- 
inary axis is considered to construct p(y), i.e. 
p{y) = Jt e dx p(x,y), where p(x,y) is the density 
of the complex eigenvalues x + iy. 

The density p(y) is used to estimate a value 
for the chiral condensate by applying the Banks- 
Casher relation (which originally was derived for 
real eigenvalues appearing in pairs of opposite 
sign). We further employed the standard defi- 
nition of the Green function j^] by inverting the 
fermionic matrix with a noisy source and by com- 
puting its eigenvalues exactly, respectively, to get 
the condensate. Thus the chiral condensate for 
two-color QCD with finite chemical potential was 
extracted by three methods. The preliminary re- 
sults for < tptp > and its variance are shown in 
the upper right plot of Fig. |[ 

In the lower left plot of Fig. || we show the den- 
sities of the small eigenvalues of the "Gor'kov op- 
erator" , averaged over 50 configurations. Again a 
band e = 0.1 is taken to construct p(y) from the 
complex eigenvalues. This density is employed in 
the spirit of the Banks-Casher relation to obtain 
a value for the diquark condensate. It compares 
well with the results of Ref. j| from an inversion 
of the K matrix (Q) with m = 0.05 and f3 — 1.5, 
see lower right plot of Fig. ||. 

2. Quantum chaos 

The fluctuation properties of the eigenvalues in 
the bulk of the spectrum have also attracted at- 
tention. It was shown in Ref. || for Hermitian 
Dirac operators that on the scale of the mean 
level spacing they are described by RMT. For 
example, the nearest-neighbor spacing distribu- 
tion P(s), i.e. the distribution of spacings s be- 
tween adjacent eigenvalues on the unfolded scale, 
agrees with the Wigner surmise of RMT. Ac- 
cording to the Bohigas-Giannoni-Schmit conjec- 
ture 0, quantum systems whose classical coun- 
terparts are chaotic have a nearest-neighbor spac- 
ing distribution given by RMT whereas systems 



whose classical counterparts are integrable obey 
a Poisson distribution, Pp(s) = e~ s . Therefore, 
the specific form of P(s) is often taken as a cri- 
terion for the presence or absence of "quantum 
chaos" . 

For p > 0, the Dirac operator loses its Her- 
miticity properties so that its eigenvalues become 
complex. We apply a two-dimensional unfold- 
ing procedure [|| to separate the average eigen- 
value density from the fluctuations and construct 
the nearest-neighbor spacing distribution, P(s), 
of adjacent eigenvalues in the complex plane. 
The data are then compared to analytical predic- 
tions of the Ginibre ensemble jt| of non-Hermitian 
RMT, which describes the situation where the 
real and imaginary parts of the strongly corre- 
lated eigenvalues have approximately the same 
average magnitude. For this case P{s) is given 
in Ref. p0| . For uncorrelated eigenvalues in the 
complex plane, the Poisson distribution becomes 
Pp( s ) = (vr/2) s exp(-7rs 2 /4) fuj. This should 
not be confused with the Wigner distribution for 
a Hermitian operator |(| . 

Our results for P(s) of the Dirac operator are 
presented in the upper row of Fig. ^. As a func- 
tion of p, we expect to find a transition from 
Wigner to Ginibre behavior in P(s) ||. For 
the symplectic ensemble of color-SU(2) with stag- 
gered fermions, the Wigner and Ginibre distribu- 
tions are very close to each other and thus hard 
to distinguish. They are reproduced by our pre- 
liminary data for p = and p = 0.4, respectively. 
Even in the deconfined phase, where the effect 
of the chemical potential might order the system, 
the gauge fields retain a considerable degree of 
randomness, which apparently gives rise to quan- 
tum chaos. 

For p > 1.0, the lattice results for P(s) deviate 
substantially from the Ginibre distribution and 
could be interpreted as Poisson behavior, corre- 
sponding to uncorrelated eigenvalues. A plausi- 
ble explanation of the transition to Poisson be- 
havior is provided by the following two observa- 
tions. First, for large p the terms containing in 
Eq. (Q) dominate the Dirac matrix giving rise to 
uncorrelated eigenvalues. Second, for large p the 
fermion density on the finite lattice reaches satu- 
ration due to the limited box size and the Pauli 
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Figure 3. Nearest-neighbor spacing distribution of the Dirac matrix (upper plots) and of the "Gor'kov 
operator" (lower plots) for two-color QCD with varying \i. The curves are the Wigner distribution Pw — 
262144/(7297r 3 )s 4 exp(-64/(97r)s 2 ), the Ginibre distribution of Rcf. Q and the Poisson distribution 
Pp(s) = (n/2)s exp(— 7rs 2 /4), respectively. 



exclusion principle. 

The results for P(s) of the "Gor'kov operator" 
are plotted in the lower row of Fig. || Since the 
eigenvalues are complex also for fi = 0, we find 
approximately Ginibre behavior corresponding to 
non-Hcrmitian RMT. For /j, = 3.0, P(s) takes a 
Poisson distribution associated with uncorrelated 
complex eigenvalues. 
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